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Two portions of a synthetic hydroxyapatite (HA), Ca>OH(P04);j, fully characterized by x-ray, 
infrared, petrographic, and chemical analyses, were heated at 1,000 °C in air and steam atmospheres, 
respectively. Solubility isotherms for these two samples in the system Ca(OH)2-H3P04-H20 were 
determined in the pH range 5 to 7 by equilibrating the solids with dilute H3PO4 solutions. Both sam- 
ples of HA dissolved stoichiometrically. The activity products (Ca ++ ) 5 (OH _ )(POf ) 3 and their standard 
errors — obtained by a least squares adjustment of the measurements (Ca and P concentrations and 
pH of the saturated solutions) subject to the conditions of electroneutrality, constancy of the activity 
product, and stoichiometric dissolution — were 3.7. } ±0.5 X 10 -58 for the steam-heated HA and 
2.5i±0.4x 10 -55 for the air-heated HA. Allowance was made in the calculations for the presence of 
the ion pairs [CaHPOd and [CaH 2 P04| + . The higher solubility product for the air-heated HA is as- 
cribed either to a change in the heat of formation brought about by partial dehydration or to a state 
of fine subdivision resulting from a disproportionation reaction. The solubility product constants were 
used to calculate the points of intersection (i.e., singular points) of the two HA solubility isotherms with 
the isotherms of CaHP0 4 ' 2H 2 and CaHP0 4 ; it was found that the pH's of the singular points for 
the air-heated HA were a full unit higher than those of the steam-heated preparation. Conditions 
are described for the precipitation of HA crystals suitable for solubility measurements. 
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1 . Introduction 

Hydroxyapatite (HA), Ca50H(P0 4 ) 3 , has wide 
biological and industrial importance. Its physico- 
chemical properties, in particular solubility and de- 
tailed crystallographic structure, are directly related 
to the mineralization processes in biological systems. 
The various reports in the literature on solubility of 
HA are contradictory. Clark [l], 1 on the basis of 
precipitation and dissolution experiments, concluded 
that HA has a definite solubility product. Thermody- 
namic functions for HA (from which a solubility prod- 
uct constant can be calculated) have also been reported 
[2, 3, 4]. Other reports [5-8], however, claim that 
systems saturated with respect to HA do not con- 
form to the principle of a solubility product constant 
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based on the stoichiometry of the salt. Among the 
models that have been advanced to explain the 
apparent anomalous behavior of HA in aqueous solu- 
tions, those advocating formation of complexes on the 
surface of the solid [5, 8] have received considerable 
attention. These models implicitly assume that 
equilibrium between the solid and the liquid phases 
is not achieved with the crystalline portion of the HA, 
but instead with surface layers having quite different 
chemical composition from that of the bulk of the 
solid. Furthermore, the composition of such layers 
has been assumed to be constant over a wide range of 
calcium, phosphorus, and hydrogen concentrations 
in the liquid phase. 

At least part ol the problem in the interpretation 
of data on solubility of HA is related to the difficulty 
in preparing pure samples of this compound. So 
far, the best crystals have been obtained by hydro- 
thermal methods, but the yields are small, the equip- 
ment is relatively expensive and the product may 
contain extraneous phases such as /3-calcium pyro- 
phosphate and/or /3-tricalcium phosphate. Direct 
precipitation methods usually yield very finely divided 
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products, often without the stoichiometric composition 
of HA. A particularly serious difficulty is contamination 
by precipitation of the more acid calcium phos- 
phates, among which octacalcium phosphate (OCP), 
Ca 8 H 2 (P04)6 * 5H 2 and its decomposition products 
are the most difficult to eliminate. 

The purposes of this work were (a) to devise a sim- 
ple method for preparation of HA that would yield 
sufficiently large and pure crystals so that an un- 
ambiguous characterization would be possible, (b) 
to learn whether this compound behaves as a normal 
chemical compound with unique thermodynamic 
properties when equilibrated with aqueous solutions, 
and (c) to establish its solubility isotherm and solu- 
bility product constant at 25 °C. 

2. Experimental Methods and 
Procedures 

In the method of preparation outlined here, the 
possibility of OCP precipitation is minimized by the 
use of relatively high temperature and pH. However, 
pH values higher than 10 are avoided because they 
seem to favor the formation of very finely divided 
precipitates. Other conditions are (a) slow addition 
of reagents to reduce continuous nucleation, (b) 
high ionic strength and the presence of acetate ions 
(a Ca complexing agent) to increase the solubility of 
HA, and thus to favor formation of larger crystals, 
(c) use of a large reaction vessel to perform the pre- 
cipitation without sacrificing yields, and (d) special 
precautions to avoid contamination with carbonate 
from the atmosphere or reagents (ammonia was used 
to control the pH instead of NaOH or KOH because 
the alkali hydroxides are difficult to obtain free of 
carbonate contamination). 

2.1. Preparation of HA 

Freshly boiled distilled water and chemically 
pure salts were used to prepare 2.5 liter of 7.25 X 10~ 2 
M (NH 4 ) 2 HP0 4 solution and 2.5 liter of 13.33 X 10~ 2 M 
Ca(NO,3)2 • 4H2O solution; their pH's were adjusted 
with ammonia gas to 8.5-9. These solutions were 
stored in borosilicate glass bottles vented through 
ascarite tubes. 

The reaction vessel was a 22-1, 3-neck, borosili- 
cate glass flask. A 4-ft condenser was attached to 
one neck; transfer of reagents and materials was done 
through the other two necks. The top of the condenser 
was connected by a glass tube to a 6-1, 3-neck flask 
containing 3 liters of 10 N NaOH solution; a check 
valve permitting outward flow of gases was con- 
nected to one neck and a tube, extending into the NaOH 
solution and acting as a relief valve, entered through 
the third neck. The reactant solutions were introduced 
into the large vessel by a peristaltic pump at a rate 
of 100 ml/hr. Ammonia gas was passed through an 
ascarite tower and then bubbled through water be- 
fore entering the reaction vessel. Samples were with- 
drawn through a tube that extended to the bottom of 



the reaction vessel for petrographic examination of 
the solid and pH determination of the liquid. 

Initially, 9 1 of A^ ammonium acetate solution were 
boiled in the reaction vessel for about 12 hr to remove 
C0 2 ; then ammonia gas was passed in until the pH 
reached a value between 8.5 and 9, and the pumping 
of the phosphate and calcium solutions was started, 
maintaining vigorous boiling. The pH was maintained 
between 8.5 and 9.5 (measured at room temperature) 
by an adequate flow of NH 3 until the total volume of 
the solutions was pumped into the reaction vessel; 
then the flow of NH3 was discontinued but boiling 
was continued for about 3 hr. The solid was allowed 
to settle, the supernatant liquid was withdrawn, 
about 9 1 of boiled distilled water were introduced 
and the mixture was boiled for about 2 hr. This 
washing operation was repeated as was deemed nec- 
essary to remove the soluble salts (usually five times). 
Further washing was done by decantation, using 
boiled distilled water and short settling periods so 
as to discard as much as possible the finest fraction 
of the precipitate. Precautions to avoid contami- 
nation by atmospheric CO2 were taken during all 
these operations. The suspension was then filtered 
and the solid was dried at 150 °C for 48 hr. As ex- 
plained in section 3, a portion of this sample was 
heated at 1,000 °C in air and another portion was heated 
at the same temperature in steam at one atmosphere. 

2.2. Equilibrations 

Equilibrations were made in glass-stoppered bot- 
tles, each of which contained 5g of HA and 120 ml 
of a phosphoric acid solution, by slowly rotating end- 
over-end in a water bath at 25 ±0.1 °C. After two 
weeks, the suspension was transferred to a thermo- 
stated saturator, described in an earlier publication 
[9], in which the pH of the clear filtered solution was 
measured. Samples of this solution were taken for 
calcium and phosphorous analyses. Standardization 
of the pH-meter (with claimed relative accuracy 
of ±0.0037 pH units) was made with certified NBS 
buffer standards [10]. 

Phosphoric acid solutions used in the equilibra- 
tions were prepared from doubly crystallized 
2H3PO4 • H 2 and freshly boiled conductance water. 
The molarities of the final solutions were determined 
by chemical analyses. 



2.3. Calcium and Phosphorus Determinations 

Calcium was determined colorimetrically by the 
indirect technique described by Banerjee et al. [11]. 
Phosphorus was determined using the vanadomolyb- 
date reagent of Brabson et al. [12]. In both analyses, 
samples were compared directly with blanks and 
standards in a double-beam spectrophotometer. 

The analytical errors estimated for the measure- 
ments reported here were: pH, ± 0.009 units; calcium 
±2.5 percent; and phosphorus, ±1.5 percent of the 
amounts analyzed. 
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3. Characterization of the Solid 

The HA was characterized by petrographic, gas- 
adsorption, chemical, infrared and x-ray methods. 

Under the petrographic microscope, the crystals 
appeared as masses of birefringent, irregularly shaped 
equant particles up to 5 /xm in diameter. Present, also, 
was a minor fraction of thin needles up to 10 /xm in 
length. The mean index of refraction (white radiation, 
25 °C) was 1.636 with a standard error of ±0.002. The 
material was petrographically homogeneous, both the 
equant particles up to 5 /xm in diameter. Present, also, 
HA. The mean index is lower than the indexes usually 
given for HA prepared hydrothermally, N € = 1.644 and 
7Vo,= 1.651 [37], but this appears to be a character- 
istic of precipitated HA. The treatments at 1,000 °C 
did not change the mean indexes significantly. 

The surface area determined by the BET procedure, 2 
using nitrogen as the adsorbate, yielded a surface area 
of 6 m 2 /g. Although this area is 5 to 10 times larger 
than that expected from microscopic examination, the 
agreement is reasonable considering the uncer- 
tainties inherent to the two methods. 

Chemical analyses yielded 18.22 ±.27 percent P 
and 39.63 ±0.99 percent Ca, the indicated errors 
being calculated on the basis of the estimated standard 
errors of the analytical procedures employed. These 
values compare with the theoretical values 18.50 
percent P and 39.89 percent Ca for HA. 

The infrared spectrum displayed all the bands 
characteristic of well-crystallized HA [14]. Small 
bands in the regions of 1410 and 1640 cm -1 indicated 
slight contamination with carbonate and moisture, 
respectively. Analysis by the Warburg method yielded 
0.04 percent C0 2 . 

The x-ray diffraction pattern was sharp and charac- 
teristic of that for pure HA [13]. No extraneous peaks 
were found. 

Although the above methods to characterize the HA 
did not reveal the presence of extraneous phases, when 
a sample of the HA was heated to 1,000 °C for 24 hr, 
the infrared spectrum showed small peaks in the 
region 920 to 1,200 cm -1 . These peaks indicate the 
presence of /3-Ca 3 (P0 4 )2 [14], which presumably 
formed on heating from a contaminant, an acid cal- 
cium phosphate. 

A portion of the HA, therefore, was heated at 
1,000 °C for 6 days in room atmosphere to anneal the 
crystals, to drive off carbon dioxide and moisture, 
and to convert the less basic calcium phosphate into 
/3-Ca 3 (P04)2. Following the heat treatment, the sample 
was extracted with 2 X 10~ 3 M H 3 P0 4 to remove the 
/3-Ca3(P0 4 )2. Approximately 10 percent of the sample 
was dissolved by this treatment. 

Another portion of the HA was subjected to the 
same heat treatment but in the presence of one at- 
mosphere of water vapor, and then extracted with the 
H3PO4 solution. 



The infrared spectra of the heated samples were 
very similar, both corresponding to pure HA [14] but 
the height of the peak at 3,570 crrr 1 was considerably 
greater for the steam-heated sample, indicating either 
a more uniform environment for the OH" group in 
this sample or actually a higher content of OH~ in 
the crystals. 

Another sample of HA, referred to here as "as 
prepared" was synthesized according to the method 
described above, but it was not subjected to high 
temperature treatment. 



4. Calculation of Solubility Product 
Constants 

The solubility product constants, 

K S7 , = (Ca ++ ) 5 (OH-)(P0 4 -)3 

where parentheses denote ionic activities, were cal- 
culated by a procedure very similar to the one de- 
scribed elsewhere [9]. The values are listed in tables 1, 
2, and 5. However, it was recognized that the best 
values to be derived from such data must take into 
account more completely the propagation of errors 
and the proper weight of the measurements. A more 
refined method was devised, therefore, for these 
calculations which entails the simultaneous ad- 
justment of the weighted observables — i.e., calcium 
and phosphorus concentrations and pH's in the various 
equilibrated solutions plus the initial concentration 
of H3PO4 — and the determination of a parameter 
which yields best estimates of K sp . The main features 
of the method follow the general least-squares adjust- 
ment discussed by Deming [15] and they are similar 
to those recently presented by Wentworth et al. 
[16] for the calculation of complexing constants 
using spectrophotometric data. However, for the pres- 
ent work, two or three condition functions were used 
instead of one. A detailed account of the method 
will be given elsewhere [17]. It suffices here to give 
the three condition functions used in this investiga- 
tion, based on electroneutrality, the solubility prod- 
uct for HA, and congruency in the dissolution of the 
solid, respectively, 
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2 Measurement made by Dr. W. V. Loebenstein, NBS. 



in which M and P are the total molar concentrations 
of Ca and P in solution, respectively; (H + ) is the ac- 
tivity of the hydrogen ion; / is the molar activity 
coefficient of the subscripted species; K w and K3 
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are the dissociation constants of water and the third 
ionization constant of H3PO4, respectively; Pq is 
the initial molar concentration of H3PO4; S is the 
parameter to be calculated, its value being the best 
estimate of K sp ; the fraction 5/3 is the theoretical 
Ca/P ratio in HA; r is the total molar concentration 
of Ca or P present as the ion pairs [CaHPC^] and 
[CaH 2 P04] + [9] and it is defined, as are the quanti- 
ties L, Q u Qi, and /?, by the expressions, 
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In eqs (4r-8), K\ and K? are the first and second ioniza- 
tion constants of H3PO4, respectively, and K x and K y 
are the association constants for the ion-pairs 
[CaHP0 4 ]° and [CaH 2 P0 4 ] + . A list of the numerical 
values for the various constants used in this investi- 
gation is given in the appendix. 

The three condition functions, eqs (1 _ 3), were 
used to obtain the best estimates of K sp for the pre- 
heated hydroxyapatites (in air and steam). Condi- 
tion (3), defining stoichiometric dissolution, was 
excluded for the calculations with HA "as prepared" 
because phosphate was adsorbed in some of these 
equilibrations, as explained in section 5.1. 

Equation (4) was obtained from the simultaneous 
solution of eqs (9) and (10) which describe the equi- 
libria involved in the formation of ion pairs. Calling 
X the concentration of [CaHP0 4 ]° and Y that of 
[CaH 2 P0 4 ] + , 



X = 



K,[M-r][P-T]fc, 



Y= 



K Q [M-r][P-T]fc a + + W 
QM 



(9) 



(10) 



in which f x and f u are the activity coefficients for the 
species [CaHP0 4 ]° and [CaH 2 P0 4 ] + , respectively. 

The ionic strength, jjl, was defined in these calcu- 
lations by 



^ = 2[M-t] + 1/2 



m + K* +RT ^-rm 



/h+ (H + )/o, 
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(H+ 



Kit 



-+ 



9K:< 



H2PO4 



./„ 



PO4 



(H + )/ f 



po 4 



(11) 
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The values for the ionic activity coefficients were 
calculated by the use of the Debye-Hiickel theory 
using the distances of the closet approach given in the 
appendix. The iterative procedure, starting with uni- 
tary activity coefficients, reported in previous publi- 
cations [9, 18] was used for these calculations. The 
activity coefficients f x and /h 3 po 4 were taken to be 
unity, and it was assumed that/y=/ H2P0 -. 

5. Results 
5.1. Solubility of HA 

The experimental results for the steam-heated and 
air-heated preparations are given in tables 1 and 2, 
respectively, along with the values of K sp calculated 
from these results using the activity product that ap- 
pears as the leading term in eq (2). The spread in K sp 
in each table is acceptably low even though the aver- 
age K sp for the steam-heated material, table 1, is three 
orders of magnitude less than that for the air-heated 
material, table 2. The other tests that can be applied 
to these data — apparent electroneutrality unbalance, 
stoichiometry of the dissolution reaction, and smooth- 
ness of the plots of Ca or P versus pH — also show 
deviations which are within the anticipated experi- 
mental errors, thus vouching for the validity of the 
results. In tables 3 and 4 are given the adjusted values 
for the three observables, Ca, P, and pH, and the 
adjusted values of K sp , electroneutrality and the 
ratio Ca/P of the dissolution reactions. It is clear 
that the adjustments in the Ca, P, and pH values are 
within their anticipated experimental errors. The fact 
that these adjustments are small is further confirmation 
of the internal consistency of the data. We believe that 
the best estimate of the two K sp s obtained by the 
adjustment procedure (tables 3 and 4), which weights 
the experimental measurements according to their 
errors, constitute better values for the solubility prod- 
ucts than the simple unweighted mean values, given 
in tables 1 and 2. 

Table 1. Solubility of steam-heated HA, unadjusted quantities 



Initial 


Composition of saturated solution 






[H,P0 4 ] 




Ca/P 


K sp X 1058 








MxlO 3 


pH 


[P] 
A/xlO 3 


[Ca] 
M x 10 3 


Dissol. Reac. 




4.17 


5.04, 


6.19 


3.16 


1.56 


3.7 


2.00 


5.40b 


2.94 


1.49 


1.58 


5.I4 


1.00 


5.69 9 


1.50 


0.79 


1.59 


5.2, 


1.00 


5.74 7 


1.48 


.78 


1.64 


9.7, 


0.30 


6.22, 


0.46 


.25 


1.50 


2.9 9 




Avg. 1.57±0.02 a 


5.3 5 ±2.6. 2 a 



a Standard errors in the means. 
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Table 2. Solubility of air-heated HA, unadjusted quantities 



Initial 


Composition of saturated solution 






(IIJ'o,] 






Ca/P 


K,pX 10" 








M x 10 :i 


pH 


[P] 
MxW 


[Ca] 
MX 10 :t 


Dissol. Reac. 




2.04 


5.76« 


3.02 


1.62 


1.65 


2.5 7 


1.98 


5.79,, 


3.00 


1.57 


1.53 


3.1 :( 


0.99 


6.11« 


1.39 


0.76 


1.90 


2.3,; 


1.01 


6.12,, 


1.46 


.83 


1 .85 


4.1., 


1.01 


6.13o 


1.51 


.84 


1.70 


5.6, 


0.50 


6.40,, 


0.75 


.49 


1.95 


3.8, 


.10 


7.07, 


.19 


.15 


L.60 


3.6„ 


.05 


7.36, 


.11 


.09 


1.62 


2.2 7 




Avg. 1.72±0.02 a 


3.4,±0.1 a 



;| Standard errors in the means. 



Table 


3. Solubility of 


steam-heated HA, adjusted quantities 




Composition of saturated 








Apparent 




Initial 


>< .hit ion 


Ca/p 




Ionic 


charge 


K* 






Dissol. 


,.a 


strength 


imbalance 


x 10™ 


Lri ;t r'U,J 








M x 10* 


pH 


LP] 
M x 10 :! 


[Ca] 
M x 10 ! 


Reac. 




ix x 10 :i 


M X 10 ,; 




4.28 


5.05 7 


6.13 


3.10 


1.68 


0.03 H 


9.10 


0.0 


3.7, 


2.00 


5.38., 


2.89 


1.47 


1.67 


.02, 


4.40 


-.1 


3.9, 


1.02 


5.69- 


1.48 


0.77 


1.67 


.01 7 


2.32 


.0 


3.9, 


1.00 


5.72,, 


1.45 


.75 


1.67 


.01h 


2.28 


-.3 


5.3, 


0.31 


6.23 ti 


0.46 


.25 


1.68 


.01, 


0.80 


.0 


3.8, 




S = 3.7 A ±0A H U 



a Fraction of Ca present as ion pairs. 

11 Standard error from least-squares adjustment. 



Table 4. 


Solubility of air-heated HA, 


adjusted quantities 




Composition of saturated 








Apparent 




Initial 




solution 


Ca/P 




Ionic 


charge 


K*p 








Dissol. 


r a 


strength 


imbalance 


x 10 58 


LHirUjJ 








M x 10 : ' 


pH 


[P] 
M x 10 :i 


[Ca] 
M x 10< 


Reac. 




ix X 10-' 


M x 10" 




2.06 


5.77, 


3.02 


1.59 


1.67 


0.03, 


4.75 


0.0 


2.5, 


2.02 


5.78 5 


2.94 


1.55 


1.68 


.03, 


4.64 


-.1 


2.5 7 


0.96 


6.11;, 


1.42 


0.78 


1.69 


.02 7 


2.39 


.0 


2.5, 


.99 


6.10, 


1.47 


.81 


1.68 


.02* 


2.48 


-.2 


2.8„ 


1.00 


6.10« 


1.49 


.82 


1.67 


.02s 


2.52 


-.5 


3.2,, 


0.50 


6.40,, 


0.77 


.44 


1.68 


.02, 


1.42 


.0 


2.6„ 


.11 


7.08,, 


.19 


.14 


1.68 


.01 7 


0.48 


.0 


2.6, 


.05 


7.39 :t 


.11 


.09 


1.67 


.0L 


0.33 


.0 


2.5, 




S = 2.5,±0.3 7 b 



'' Fraction of Ca present as ion pairs. 

h Standard error from least-squares adjustment. 



The large difference between the two K sp values 
(tables 3 and 4) probably relates to structural changes 
induced by the heat treatments, as discussed later. 

Also listed in tables 3 and 4 are the values for 
ionic strength and the fraction, r, of Ca present in the 
form of ion pairs given by r/M. The gradual increase 
in r with decreasing pH relates to the increase in Ca 
and P concentrations, and is mostly the result of an 
increase in the concentration of the ion-pair 
[CaH 2 P0 4 ] + . 

Results of the solubility experiments for the "as 
prepared" HA are shown in table 5. The large vari- 
ation in the K sp values in table 5 is believed to relate 
to the presence of other calcium phosphates less basic 



than HA; evidence for this type of contamination is 
presented in section 3. The data also suggest the 
presence of an anion foreign to the system 
Ca(OH)2-H 3 P04-H 2 as evidenced by an apparent 
positive electroneutrality unbalance which is far be- 
yond experimental error for the three treatments 
with pH's above 6. This may have been caused by 
carbonate since the infrared spectrum showed small 
CO;7 peaks. 



Table 5. Solubility of "as prepared" HA, unadjusted quantities 





Final solution composition 




Initial 










[H,PO,] 








K<„x 10 58 


M x 10 ' 


pH 


IP] 
4/X10' 


[Ca] 

Mx 10 ' 




7.82 


4.74 H 


10.17 


5.11 


0.8, 


3.05 


5.26s 


3.43 


1.67 


1.5, 


0.49 


6.60., 


0.23 


0.16 


L1.9< 


.20 


7.05, 


.06 


.10 


10.9, 


.05 


7.45« 


.03 


.11 


251.2 



Interpretation of the data in table 5 is further 
complicated by sorption phenomena; for example the 
phosphorus concentrations in the solutions for the 
last three entries of table 5 are less than the initial 
H3PO4 concentrations. A systematic study of this 
type of surface phenomenon with finely divided HA 
has been reported elsewhere [20]. Although some of 
the K S p values in table 5 are in the range of those 
given in table 1 for the steam-heated HA, the three 
at the higher pH's are abnormally large. Furthermore, 
the calculated adjustments of the measurements, not 
shown here, were so large as to invalidate the least 
squares procedure. As a result, little weight is given 
to the set of data in table 5. They give some confirma- 
tion in the low pH range of the results obtained for 
the steam-heated HA and they illustrate how the solu- 
tions can be made supersaturated in the high pH range 
by the presence of calcium phosphates less basic 
than HA. 

The phosphorus concentrations of the solutions 
saturated with respect to the preheated HA's, tables 
1 and 2, are plotted against their corresponding pH's 
in figure 1. Similar plots for solutions saturated with 
respect to CaHP0 4 ■ 2H>0 (DCPD) [17, 18] and 
CaHP0 4 (DCPA) [19] at 25 °C are also shown in 
the same figure. The HA curves were extrapolated 
(dotted lines) to meet the intersection points A, B, C, 
and D in figure 1. These four points correspond to 
singular points between the preheated HA's and DCPD 
or DCPA. The chemical compositions of the solutions at 
the singular points were calculated by solving eq (13) 
for (H+) through a method of successive approximations, 



2(H»*k dc k | (H+) | rk,k dc 



A 



Ca + + 



/m 



R 



K w 



(H+YK (H + )/ 



-=0 



(13) 
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pH 

FIGURE 1. Solubility isotherms of steam-heated HA, full circles, and air-heated HA, 

open circles, at 25 °C. 

The two lighter curves are the isotherms of CaHP0 4 • 2H 2 [17], DCPD, and CaHP0 4 [19], DCPA. Compositions 
of the calculated singular points A, B, C, and D are given in table 6. 



in which K DC represents the solubility product constant 
for either DCPD or DCPA, as the case may be, and 
K is defined by 



K = 



K« 



K l K * K D( 



(14) 



The value for (H + ) obtained from eq (13) was then 
used to obtain the numerical value of t since, for 
saturation with respect to either DCPD or DCPA, 



K X K 



DC 



l-R • 



(15) 



The concentrations of calcium and phosphorus are 
obtained from eqs (16) and (17), respectively. 



M= wrK K K 



./; 



Ca + 



(E + fK 



(16) 



(17) 



The chemical compositions of the solutions at the four 
singular points, calculated in the foregoing fashion, 
are given in table 6. Considerably higher concentra- 
tions of Ca and P (and lower pH's) are found at the 
singular points of the preheated HA's with DCPD than 
at the corresponding singular points with DCPA; the 
same pattern is observed for the fractions of Ca present 
as ion pairs (r in the table). These results are a conse- 
quence of the higher solubility of DCPD [17, 18] 
relative to the anhydrous salt at 25 °C [19]. The most 
striking differences in the compositions at the singular 
points, however, are brought about by the solubility 
of the HA itself. Thus, a full pH unit difference exists 
between the singular points of air-heated and steam- 
heated HA with either of the two dicalcium phosphates; 
changes in calcium and phosphorus concentrations 
by factors of about 3.5 and 4, respectively, parallel 
this pH difference. If, as it is apparent from the present 
results, the solubility of HA depends upon the previous 
history of the sample, it is nearly impossible to predict 
the composition of singular points between a precipi- 
tated HA and other calcium phosphates without fully 
characterizing the solubility properties of the HA 
involved. 
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Table 6. Solution compositions at singular points 





Solid phases 


pH 


[P] 
MxlO 3 


[Ca] 
MxlO 3 


r A 


DCPD + HA (air) 


5.31 5 
4.41o 
5.70, 
4.79, 


8.64 
30.21 

3.53 
11.54 


4.43 

15.09 

1.85 

5.80 


0.05 4 


DCPD + HA (steam) 


,10b 


DCPA + HA (air) 

DCPA + HA (steam) 


.03b 
.05* 







a Fraction of Ca present as ion pairs. 



Although the weight loss in our sample was not 
measured, slight weight losses have been reported 
[22] in the range 1,000 to 1,100 °C for synthetic HA 
samples heated gradually from 25 to 1,200 °C; these 
losses occurred without alteration in the x-ray 
patterns. Also, it has been reported [23] that at 1,200 °C 
in air, disproportionation of HA into tricalcium phos- 
phate and tetracalcium phosphate occurs according 
to the reaction, 



2Ca :) OH(PO|): } -> 2 a-Ca 3 (P0 4 ) 2 



+ Ca 4 0(P0 4 )> + H>0. 



(18) 



The x-ray powder diffraction pattern of the air-heated 
sample was the same as that of HA, but it is impossible 
to rule out a partial disproportionation. 

In view of the foregoing considerations, there are 
two explanations for the higher solubility of the air- 
heated HA: (1) Reversal of the disproportionation 
reaction [18] when the sample is returned to an aqueous 
system might yield very finely divided HA with an 
unusually high solubility because of surface energy 
contribution, or (2) loss of OH" ions by dehydration 
during heating of the solid, without destruction of the 
HA lattice, might yield crystals with higher free energy. 
Evidently, the process responsible for the increased 
solubility is not readily reversible, since otherwise 
the values for the K sp of the two preparations would 
have been the same. Therefore, the equilibria attained 
with the air-heated HA must have been metastable 
relative to the more stable form of HA. This would 
require that the rate of dissolution of the air-heated 
HA was higher than the rate of precipitation of the 
more stable salt. 

The K sp values for the two kinds of preheated 
apatites correspond to a difference of about 16 kj/mol 
(3.9 kcal/mol) in the standard free energies of forma- 
tion. If it were assumed that the standard heat of 
formation is the same for both compounds, the ob- 
served difference in free energy would correspond to 
about 59 J/mol K (14 cal/mol K) less entropy in the 
air-heated apatite than in the steam-heated salt. 
This difference appears too large to be associated 
with a change in entropy resulting from vacancies 
or ordering. For this reason it is believed that the 
difference in solubility is mostly related to a difference 
in their heat contents. 



6. Discussion 

6.1. Solubility of HA 

The solubility data given in tables 1 and 2 for the 
two preheated HA samples are consistent with 
stoichiometric dissolution represented by the reaction 
Ca 5 OH(P04). } ^5Ca + + -hOH-h3POf. It is believed 
that saturation conditions were obtained in the two- 
week equilibrations as evidenced by the smooth curves 
obtained in the plots of calcium concentration versus 
pH (fig. 1), calcium concentration versus phosphorus 
concentration (not shown in this paper), and a rea- 
sonable constancy in the ionic activity product for 
HA. The success of the adjustment procedure de- 
scribed in section 4 is additional evidence that the 
parameters S, tables 3 and 4, are valid solubility 
product constants for the preheated hydroxyapatites. 

The higher solubility product constant obtained for 
the HA heated in air, relative to the product obtained 
for the sample heated in steam, is most probably due 
to structural changes induced by the thermal treat- 
ment. Such changes might relate in part to the statis- 
tically disordered structure of HA described by Kay 
et al. [21], originating from small displacements of the 
OH~ groups from the mirror planes in the space group 
P6 ; j/m. If this were purely an entropy effect, such 
disorder would tend to stabilize the lattice; the steam- 
heated sample, therefore, would have to be the dis- 
ordered one. However, it is difficult to see why heating 
HA in air would induce more order in the OH~ groups 
than would be obtained by heating in steam. A more 
likely explanation is that heating in air results in 
loss of OH~ groups by partial dehydration of the HA. 
Indeed, the infrared spectra tend to substantiate this 
alternative. Whereas the half width of the peak at 
3,570 cm -1 corresponding to the OH" stretching 
mode, was essentially the same for both preheated 
hydroxyapatites, its height for the steam-heated HA 
was almost twice that of the air-heated HA. 

6.2. Singular Points 

A system of two solid phases in equilibrium wit 1, 
an aqueous phase in the ternary system Ca(OH) 2 - 
H.SPO4-H2O is characterized by zero degrees of 
freedom (under conditions of constant temperature 
and pressure). Such singular points are important 
in solubility diagrams because they correspond to 
points where there is a change in the relative stability 
of the two solid phases. Thus, with reference to table 
6, in saturated solutions with pH's below 4.33, DCPD 
is more stable than the steam-heated HA; in solutions 
with pH's above 4.33, steam-heated HA is more 
stable than DCPD. The greater solubility of the 
air-heated HA has the effect of increasing the pH of 
the singular points with the dicalcium phosphates 
by about one pH unit relative to those with steam- 
heated HA. If an apatite comparable to the air-heated 
HA were formed in oral environments, its greater 
solubility would extend the stability ranges of the 
two dicalcium phosphates to conditions which occur 
in saliva in vivo. 
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In the concentration ranges covered by the present 
experiments, the degree of ion-pair formation varied 
from about 1.4 percent (at the highest pH) to about 
3.4 percent (at the lowest pH) of the total calcium 
concentration. In view of the uncertainty, both theo- 
retical and experimental in these calculations for 
nonsymmetrical electrolytes, the formation of ion- 
pairs in the experimental segments of the solubility 
isotherms for HA is not considered significant. How- 
ever, in more acid solutions, such as those of the 
calculated singular points, the ion-pair formation 
does become significant. 

6.3. Stoichiometry and Potential Diagrams 

The stoichiometry of a solid calcium phosphate 
equilibrated with aqueous solutions is easily revealed 
by the use of "potential diagrams" [24, 25, 26]. These 
are solubility diagrams in which the coordinates are 
functions of the chemical potentials of the components, 
Ca(OH) 2 and H3PO4, in solution. For example, if the 
quantity —log (Ca ++ )(OH~) 2 is plotted against —log 
(H + ) 3 (POj), the saturation condition is defined by a 
straight line and the absolute value of its slope cor- 
responds to the reciprocal of the Ca/P ratio in the 
solid. This kind of diagram is illustrated in figure 2 
with calculations based on the data for the two pre- 
heated apatites (tables 1 and 2). Similar calculations 
were made with data from other investigators and 
these are plotted in the same diagram. The slopes 
and their standard errors, of the regression lines for 
the air-heated and steam-heated apatites shown in 
figure 2 correspond to apparent Ca/P ratios of 1.69 



and 1.68 with standard errors of 0.02 and 0.07, respec- 
tively. These ratios agree with the theoretical ratio 
for HA, 1.67, and the ratio, 1.68 obtained by chemical 
analysis of the untreated solid. In view of these facts, 
the relatively low value obtained for the Ca/P ratio 
in the dissolution reaction, 1.57, of the steam-heated 
preparation, table 1, must be ascribed to experimental 
errors. It is to be noted that the Ca/P ratio calculated 
from dissolution involves a difference between two 
experimental measurements. 

6.4. Comparison With Reported Data 

Clark [1] reported a value of 57.75 as his best 
estimate of pK sp (the negative logarithm of K sp ) of 
HA at 25 °C. (The actual value he reported, 115.5, is 
based on the formula Caio(OH) 2 (P04)« rather than the 
formula Ca-,OH(P04)3 used here.) His systems con- 
sisted of HA precipitated by mixing dilute solutions 
of Ca(OH) 2 and H.,P0 4 at 90 °C; he kept the resulting 
suspensions at that temperature for 120 hr and 
reequilibrated them at 25 °C for 90 hr. The K sp cor- 
responding to his pK sp value of 57.75 is 1.78 XlO -58 , 
which is in good agreement with the value 3.73 X 10~ 58 
reported here for the steam-heated HA. The individual 
values reported by Clark range from 0.53 X 10 -58 to 
7.59 XlO" 58 . Even though Clark's data, as plotted in 
the chemical potential diagram, figure 2, have some 
scatter, it is apparent that they approximate closely 
the same line as our data for the steam-heated HA. 
The same author observed that equilibrations made 
at 25 °C with HA precipitated at 40 °C resulted in 
large apparent supersaturations with respect to the 




•log (H + ) (POi 



et al 



FIGURE 2. Potential plot, calcium hydroxide versus phosphoric acid. 

• steam-heated HA; O, air-heated HA; O, data from Clark [1]; A , data from Levinskas and Neuman [31]; □, data from Rootare 
[5] at 25 °C; ■. data from Rootare et al. [5] at 40 °C. 



780 



values for K sp given above. In light of present knowl- 
edge, one would expect that those precipitates con- 
tained some of the more acid calcium phosphate which 
are more soluble than HA. Thus, the more reliable 
of Clark's results agree with those reported here for 
steam-heated HA. We consider it significant that such 
good agreement was obtained between materials 
prepared by precipitation from aqueous systems and 
a material subjected to annealing at 1,000 °C in a 
steam atmosphere. This suggests that the K sp value 
obtained for the steam-heated material is characteristic 
of any well-crystallized HA. 

A solubility product constant for HA was calcu- 
lated from the following thermodynamic quantities: 
A#/ = - 13,439 kj/mol (-3,212 kcal/mol) reported by 
Smirnova et al. [2] for HA which had also been heated 
at 1,000 °C in the presence of steam; S°= 780.7 J/mol K 
(186.6 cal/mol K) reported by Egan and Wakefield 
[4] for another HA treated in a similar way, and 
standard entropies of formation of the elements and 
standard free energies of formation of ions reported 
in NBS Publications [27, 28]. The product thus ob- 
tained, based on the formula CasOH^O^, is of the 
order 10" 57 . 

This result, too, is considered to be in good agree- 
ment with the order of magnitude given here, 10~ 58 , 
for the steam-heated apatite, considering the un- 
certainties in some of the thermodynamic quantities 
used and the great effect that small variations in these 
quantities have on the K sp calculated in this fashion. 

The foregoing results are not in agreement with 
those of Levinskas [29], Levinskas and Neuman [6], 
and Neuman and Neuman [7] who concluded that 
synthetic preparations of HA equilibrated in aqueous 
solutions do not obey the solubility-product constant 
principle. One of the hydroxyapatites (L-apatite) most 
extensively studied by these authors was a commercial 
preparation of fine particle size, 68 m 2 /g. The large 
surface areas of their materials could introduce dif- 
ficulties in the interpretation of solubility studies. 
However neither theoretically [30] nor experimentally 
[20] do sorption phenomena invalidate the application 
of the solubility product constant to systems equili- 
brated with HA. It should be pointed out that the 
sample of HA (L-apatite) used by Neuman and co- 
workers had a refractive index of 1.605 [5] which is 
much too low for a well-crystallized HA; in addition, 
their preparation was contaminated with calcite, 
contained 0.28 percent sulfate and 0.8 percent C0 2 
[5]. No history is given for the synthesis of L-apatite, 
but the analytical results clearly show the heteroge- 
neity of the sample. Under these circumstances, 
reactions other than simple dissolution and precipita- 
tion of HA are to be expected, and their involved 
kinetics is a burden for sound interpretations of 
solubility studies. The adoption of the ion-activity 
product (Ca+^KHPOj) to characterize HA systems, 
as suggested by Neuman and co-workers [7, 29] 
cannot be justified on theoretical bases [30]. Our 
recalculation of their data [31] plotted in figure 2 
shows that the points are widely scattered, do not 



define a single line, and cannot be used to derive a 
reliable K sp for a solid of any Ca/P ratio. 

Rootare, Deitz, and Carpenter [5] reported results 
on solubility studies done with several commercial 
and laboratory-prepared samples. They observed wide 
variations in the calculated ion-activity products for 
HA; these were explained by postulating formation 
of a surface complex according to the reactions, 

Ca 1 o(P0 4 ) <J (OH) 2 + 6H 2 0^4Ca 2 (HP0 4 )(OH) 2 ] 

+ 2Ca + ++2HP0 4 = (19) 



Ca 2 (HP0 4 )(OH) 2 *± 2Ca + + + HPOj + 20H 



(20) 



It is true, as they claim, that the ionic product cal- 
culated according to eq (20) showed less variability 
(pK values from 24.6 to 26.8) than those calculated for 
HA on the basis of the formula unit Caio(OH} 2 (P0 4 )6 
(pK from 108 to 120). However, the mechanism postu- 
lated by Rootare et al. is open to the following criti- 
cisms: (a) Although the composition of the HA surface 
does not have to conform necessarily to the stoichiom- 
etry of the bulk solid, there is no a priori reason to 
postulate that certain surface compositions will 
prevent equilibrium with the crystalline HA lattice 
[30]; (b) It is highly improbable that the surface of 
HA equilibrated in aqueous solutions has a definite 
composition independent of the composition of the 
aqueous medium [20]; (c) If a definite composition 
is ascribed to the "surface complex" and the solubility 
product principle is applied according to its stoichiom- 
etry, it is mandatory that such species be considered 
as another solid phase (i.e., has constant molar free 
energy and stoichiometry [30]); then, in the light of 
eqs (19) and (20), the system would have zero degrees 
of freedom (in the system Ca(OH) 2 -H 3 P0 4 -H 2 0) under 
fixed temperature and pressure [30] so that the solu- 
tion composition could not vary as observed by them; 

(d) It can be easily shown that the equilibrium con- 
ditions [19] and [20], as given by Rootare et al., also 
imply equilibrium with respect to the HA solid phase; 

(e) Of the 12 samples used by Rootare et al., only four 
were petrographically homogeneous [5]; the rest were 
contaminated with calcite varying from traces to about 
1 percent, with sulfate varying from 0.006 to 0.29 
percent and, with the exception of two, all the samples 
contained carbonate in significant amounts; only the 
four samples prepared in the laboratory showed a 
refractive index approaching the one expected for 
well-crystallized HA; one of the latter, P-20, probably 
contained less basic calcium phosphates as evidenced 
by a Ca/P ratio of 1.583 and the report by Gee [32] 
that pyrophosphate was formed upon heating it at 
600 °C for 60 hr. 

The results reported by Rootare et al. [5] form two 
distinct populations when plotted in figure 2. The 
regression line for the points at 25 °C (the point 30.20, 
19.61 excluded) has a slope of —0.778, which would 
correspond to a ratio of Ca/P= 1.28 in the solid phase. 
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The points from treatments at 40 °C appear to define 
two lines with slopes — 0.309 and — 0.211 which would 
correspond to Ca/P ratios of 3.24 and 4.74 in the solid. 
If equilibrium with respect to Ca 2 (HP0 4 ) (OH> 2 had 
been established, as claimed by Rootare et al. [5], 
lines with slope of —0.5 (Ca/P = 2) should have been 
defined by their data in the potential plot. The fact 
that their data define two or three Ca/P ratios, none of 
which correspond to their postulated surface complex, 
along with the other criticisms listed above, argue very 
strongly against the existence of the surface com- 
plex hypothesized by them. LaMer [8] reexamined the 
data of Rootare et al. [5] and came to the conclusion 
that the product, (Ca + +) (HPO=) showed the least 
variation. However, for the reasons given above, this 
is not an adequate basis for any conclusions regarding 
the nature of the equilibrating solid phase. 

In the light of the foregoing results and discussion 
we conclude that our samples of preheated HA did 
equilibrate according to the stoichiometry of HA and 
according to the usual solubility product principles, 
that the steam-heated HA yielded a solubility constant 
usable as a frame of reference for well-crystallized 
HA, and that most of the difficulties reported in the 
literature on this matter relate to contamination of 
the solid and lack of equilibrium between the solution 
and the solid. 



7. Appendix 

Ionic activity coefficients were calculated from the 
extension of the Debye-Huckel limiting law, 



log/i : 



AZ'fVjl 
l + Ba,V/I 



(21) 



The values used [33] for the temperature-dependent 
constants A and B were 0.5092 and 0.3286 X 10 8 , 
respectively; those for the parameters a% [34] were 
Ca + + , 6X 10- 8 cm; H 2 P0 4 -, HP0 4 =, and POf , 4X 10~ 8 
cm; H + , 9X10 8 cm; and OH~, 3.5 X 10 8 cm. 

The values used for the three ionization constants 
of H3PO4 were K u 7.10 8 X10" 3 [35]; K 2 , 6.31 8 X 10~ 8 
[36]; Ks, 4.73 X 10 13 [37]; and for K w , 1.012 XlO" 14 
[38]. The association constants used in this paper 
are defined by 



K x 



[CaHPQ 4 ]° 
(Ca++)(HP0 4 = 



K y = 



(CaH 2 PO+) 
(Ca++)(H 2 P0 4 ) 



with values [17] K x , 244.0; and K y , 8.20. 

The points in the potential plot (fig. 2) corresponding 
to the data of Rootare et al. [5], at 40 °C were cal- 
culated using the following constants: K u 5.99 6 X 10~ 3 
[35]; K 2 , 6.60 7 X10" 8 [36]; K 3 , 7.58 6 X 10~ 13 [39]; 
K w , 2.936 XlO" 14 [38]; A, eq (21), 0.524! [33]; and 
B, 0.3318 XlO 8 [33]. 
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